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Abstract 

The charm quark contribution to the first moment of gi(x,Q 2 ) is calculated 
using a heavy mass expansion of the divergence of the singlet axial current. It is 
shown to be small. 
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The size of a possible intrinsic charm contribution in the proton has been the topic 
of intensive discussions [p], |2|, [3|, |J for many years. It is therefore a natural question 
to investigate the polarized intrinsic charm distribution in the nucleon ||. Recently 
one of us and collaborators argued that earlier treatments of the polarized charm 
contribution to the rj and rf [|7|, |8| were incorrect. In this contribution we extent and 
adopt that analysis to the nucleon. More precisely we shall focus on the intrinsic 
charm contribution to the first moment of the spin structure function gi(x,Q 2 ). This 



is known to be intimately related to the gluonic axial anomaly |j [K| [11], [T2|| . It may 
be expressed as forward limit of (t), the form factor in the proton matrix element 
of the singlet axial current 

(N(p 2 ,X 2 )\j^(0)\N( Pl ,X 1 )) 

= ufr\ P2 ) (Gf(t) %l5 - G$>(t)qsy B ) u^\ Pl ), (1) 

where q = p 2 — pi and t = q 2 . The singlet pseudoscalar form factor does not acquire 
a Goldstone pole at t = 0, even in the chiral limit, contrary to the matrix elements of 
the octet currents. In this limit, there exist eight massless pseudoscalar mesons serving 
as Goldstone bosons. However, the ninth pseudoscalar, the 7/-meson, remains massive, 
due to the mixing with the QCD ghost pole. 

This fact allows to relate the forward matrix element of the axial current to the 
(slightly) off-forward one of its divergence: 



\im(N(p 2 ,X 2 Wji!l(0)\N(p 1 ,X 1 )) 



2m N G A (0)u { i 2 \p 2 ) l5 u { i 1 \ Pl ), (2) 



m^v being the proton mass. The divergence of the singlet axial current in turn contains 
a normal and an anomalous piece, 

= 2* £ m q q l5 q - (^) G%G^\ (3) 

where Nf is the number of flavours. The two terms at the r.h.s. of the last equation 
are known to cancel in the limit of infinite quark mass O, O, O] . This is the so-called 



cancellation of physical and regulator fermions, related to the fact, that the anomaly 
may be regarded as a usual mass term in the infinite mass limit, up to a sign, resulting 
from the subtraction in the definition of the regularized operators. 
Consequently, one should expect, that the contribution of infinitely heavy quarks to the 
first moment of g\ is zero. This is exactly what happens in a perturbative calculation of 
the triangle anomaly graph ||12|| . One may wonder, what is the size of this correction for 



large, but finite masses and how does it compare with the purely perturbative result. 
To answer this question, one should calculate the r.h.s. of ([3D for heavy fermions. The 
leading coefficient is of the order m~ 2 , and its calculation was addressed recently by 
two groups 0] and M who came up with results differing by a factor of six. However, 
the operator fabcG a ^ v G b va G c a ^ appearing in both treatments, does not satisfy some basic 
properties, such that both calculations seem to be flawed. 
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i) It is not a divergence of a local operator, therefore it is not clear that its forward 
matrix element @ will vanish. 

iii) It makes no contact with the calculation of the triangle diagram in momentum 



space fllj, [12] being essentially non-abelian. 



The recent contribution || corrected this result and arrived at the expression: 

• c _ a. 



5 % = t^-^r, (4) 



where 

R, = d, (G^GT*) - 4 (D a G ua ) a G% . (5) 

[Here we use the conventions: 75 = i7°7 1 7 2 7 3 and £0123 = 1] This result is an explicit 
4-divergence and has the Abelian limit. Moreover, this result can also be obtained 
by a 1/m expansion of the triangle diagram contribution. Actually the result of 
demonstrates that in order m~ 2 the entire result (y) can be restored from the venerable 
triangle diagram. The diagrams with larger number of "legs" give only contributions to 
the non-abelian part of the result Indeed, computing the forward matrix element 
of operator (|5|) between two virtual gluon states we get the following expressions: 

(p\t^- 2 RM = -*^w^v4 • ( 6 ) 

On other hands the result of a calculation of the triangle diagram with massive fermions 
(see e.g. [[□]]) has a form: 

(P\cwc\p) = i^^V{l-f^^%4} 
2ir L Jo mi — p z x(l — x) J 



p 2 x(l — x) 

-^W"eV^ + 0(ij). (7, 



This expression coincides exactly with the result (0). In order to complete the proof it 
is enough to consider the off-forward matrix element of the operator (|4]) between two 
gluons at zero virtuality and compare the result with the expression for the triangle 
diagram for (p'\ d^c^^y^c\p) . It is easy to check that again the results coincide. 
The proton matrix element of R^ takes a form analogous to that of (0) 

(#(^2)1^(0)1^1, AO) 

= u^\p 2 ) (G^(t) 7M 75 - G*(t)q» l5 ) u^\ Pl ), (8) 

It is crucial, that because of the explicit gauge invariance of R^ the zero mass ghost 
pole does not contribute. This make an apparent difference with respect to the massless 
case, when the divergence of the gauge- dependent topological current appears and 
the ghost pole contribution does not allow to deduce the relation between the matrix 



elements of the currents starting from the relation for their divergencies [13]. In the 
case under investigation only the contribution of the massive 77 meson may appear so 
that 



lim(#(p 2) A 2 )|^(0)|#(p 1) A 1 )) 



2m N G R A ^)u^\p 2 ) lb u^\p 1 ), (9) 
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The contribution of charm to the forward matrix element can be obtained by substi- 
tuting (|l], H) into the proton matrix elements of (f|), giving in the forward limit. 

(N(p,\)\j$(0)\N(p,\)) = J^(N( P ,\)\R,(0)\N(P, A)) (10) 

In deriving this expression we used |^). Note that the first term in does not 
contribute to the forward matrix element because of its gradient form, while the con- 
tribution of the second one is rewritten, by making use of the equation of motion, as 
matrix element of the operator 

(N(p,X)\3 { S(0)\N(p,X)) = t£^ c (N(p,\)\9 Ef=uAsM^V^|iV(p,A)> 

= 12wm^ m N S lifs > (11) 

The parameter f s was determined before in calculations of the power corrections to 
the first moment of the singlet part of gi part of which is given by exactly the quark- 
gluon-quark matrix element we got. Note that within our l/m c approximation the 
c contribution to the flavour sum can be neglected. QCD-sum rule calculations gave 
fs* = | U {2) (proton) + /( 2 ) (neutron)) = 0.09Q |T§, estimates using the renormalon 
approach led to f s = ±0.02 [jnj and calculations in the instanton model of the QCD 



vacuum give a result very close to that of QCD sum rule |l5j \ $% . 



Inserting these numbers we get finally for the charm axial constant the estimate 

(0) = ffC^f « -5 • 10- 4 (12) 



with probably a 100 percent uncertainty (see e.g. |18[j). As the mass term in the 



triangle diagram is coming from the region of transverse momenta of the order m c , this 

(2) 

should be the correct scale of both a s and f s . Because this scale is not far from the 

(2) 

typical hadronic scale at which fg was estimated we can neglect evolution effects. 
Note that this contribution is of non-perturbative origin (therefore we call it intrinsic), 
so that it is sensitive to large distances, as soon as the factorization scale is larger than 
m c . If the scale is also larger than m&, one can immediately conclude that the non- 
perturbative bottom contributions is further suppressed by the factor {m c jm^) 2 ~ 0.1. 
Let us note, that the naive application of our approach to the case of strange quarks 
gives for their contribution to the first moment of g\ roughly — 5-10 -2 , which is compati- 
ble with the experimental data. The possible applicability of a heavy quark expansions 
for strange quarks in a similar problem was discussed earlier JTPJ in the case of the 



vacuum condensates of heavy quarks. That analyses was also related to the anomaly 
equation for heavy quarks, however, for the trace anomaly, rather than the axial one. 

Let us summarize: We have related the non-perturbative contribution of charm quarks 
to the nucleon spin (at scale m c ) to the singlet twist-4 coefficent appearing e.g. in 
the Ellis- Jaffe sum rule. Numerically it is found to be very small, contrary to the 
suggestion of |5|, |J. We would like to note that in a recent paper [^UJ it was shown that 



Note that here we use a convention for the e-tensor which differs by sign from that of fll5| 
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also the perturbative Ac contribution is very small. We see this as further support for 
our result. 
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